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ABSTRACT 

We have used high precision differential astrometry from the Palomar High-precision Astrometric 
Search for Exoplanet Systems (PHASES) project and radial velocity measurements covering a time- 
span of 20 years to determine the orbital parameters of the 88 Tau A system. 88 Tau is a complex 
hierarchical multiple system comprising a total of six stars; we have studied the brightest 4, consisting 
of two short-period pairs orbiting each other with an '^18-year period. We present the first orbital 
solution for one of the short-period pairs, and determine the masses of the components and distance 
to the system to the level of a few percent. In addition, our astrometric measurements allow us to 
make the first determination of the mutual inclinations of the orbits. We find that the sub-systems 
are not coplanar. 

Subject headings: techniques:interferometric-star:88 Tau 
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1. INTRODUCTION 

88 Tau (HD 29140, HR 1458, HIP 21402) is a bright 
[mv = 4.25, = 3.69 ± 0.25; Skrutskie et al. 
2006) , nearby (~ 50pc ) hierarchical sextuple stellar sys- 
tem (jTokovininI Il997f ) . The A component contains a 
pair of systems (designa ted Aa and Ab) in an 18- 
year (iBa lega et al.l Il999l ) orbit that has been r esolved 
by speckle interferometry (jMcAlister et al.l 119871 ). The 
Aa component is a known spectroscopic binary system 
iP ~ 3.57-day), with a composite spectral type of A5m 
(iCowlev et all [l969l). In previous work it had been 
noted ( Burkhart &: Couprvlll988f ) that the A system is 
likely complex, wit h possibly as many as 5 components. 
iBalega et al.l ()1999| ) noted a discrepancy between the to- 
tal estimated mass of this system based on photometry 
and spectral types, and the total mass derived from the 
visual orbit and Hipparcos parallax. In this work we have 
determined that, like the Aa component, the Ab compo- 
nent is a double-lined binary; this newly-resolved binary 
has a period of 7.89 days. Finally, there is a common- 
proper-motion companion, labeled B, located ~69 arc- 
seconds awa y from the A sy stem; it, too, is known to 
be a binary (iTokovinin fc Go rvnva 2001j) . For clarity we 
provide a schematic of this complex system in Figure [TJ 

There are several reasons why multiple stellar systems 
such as 88 Tau merit attention: first, binary orbits make 
it possible to measure accurate stellar masses and dis- 
tances, while the larger number of presumably co-eval 
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Fig. 1. — A schematic diagram of the 88 Tau system. 



stars allows one to impose the additional constraint that 
any given model must accurately match all of the stars. 
This approach has proven particularly fruitful when ap- 
plied t o another fa mous hierarchical sextuple system: 
CastordTorres fc Rib as 2002). Second, as outhned in 
ISterzik fc TokovininI (|2002). the relative orientations of 
the orbital angular momenta allow one to constrain the 
properties of the cloud from which the stars are thought 
to have formed, as well as the subsequent dynamical de- 
cay process. Despite their value, observational problems 
have limited the number of triple or higher-order systems 
with accurately measured orbits to fewer than 10. Given 
their hierarchical nature it is often the case that either 
the close system is unresolvable or the outer system has 
an impractically long orbital period. 

With the advent of long-baseline stellar interferom- 
etry, and more recently p hase-referen ced long-bascline 
interferometric astrometry (jLane fc Muterspaughll200l 
capable of 10-20 /.t-arcsecond astrometric precision be- 
tween pairs of stars with separations in the range 0.05- 
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1 arcsecond, it has become possible to resolve the 
orbital motion of several interesting multiple systems 
(|Muterspaugh et al]|2006bl la[). Here we report on astro- 
metric and radial velocity measurements of the 88 Tau A 
system, which allow us to constrain the orbits of the 3.57- 
day, 7.89-day and 18-year components with improved 
precision, and for the first time provide a relative ori- 
entation of the orbits as well as component masses. 

Astrometric measurements we re made with the 
Palomar Testbed Interferometer fColav ita et al.l |1999|) 
as part of the Palomar High-precision Astrometric 
Search for Exoplanet Sy stems (PHASES) program 
(jMuterspaugh et al.l [2006d ) . The Palomar Testbed In- 
terferometer is located on Palomar Mountain near San 
Diego, CA. It was developed by the Jet Propulsion Lab- 
oratory, California Institute of Technology for NASA as 
a testbed for interferometric techniques applicable to the 
Keck Interferometer and the Space Interferometry Mis- 
sion (SIM). It operates in the J (1.2 /zm), H (1.6 /xm), 
and K (2.2 /im) bands and combines starlight from two 
out of three available 40 cm apertures. The apertures 
form a triangle with 86 and 110 m baselines. 

2. OBSERVATIONS & MODELS 

2.1. PHASES Astrometry 

88 Tau A was successfuly observed with PTI on 29 
nights in 2004-2007 wit h the use of the phase-refe renced 
fringe-scanning mode ( Lane fc Muterspaugh|[200^ devel- 
oped for high-precision astrometry; the data were re- 
duced with the algorithms describe d therein, as well as 
with the modifications described in IMuterspaugh et aTl 
([2005). 

The obtained differential astrometry is listed in Table 
[TJ Note that the astrometry on any single night is es- 
sentially that of a single-baseline interferometer, yield- 
ing a very small error in the direction aligned with 
the baseline, but limited to the effect of Earth-rotation 
synthesis in the perpendicular direction. The median 
minor-axis formal uncertainty is 10 /iarcseconds, while 
the median major-axis uncertainty is 312 /iarcseconds. 
To properly weight the data set when doing a com- 
bined fit with previous astrometry and radial velocity 
data, we fit an orbital model to the PHASES astrome- 
try by itself, and rescaled the formal uncertainties so as 
to yield a reduced of unity; the resulting scale fac- 
tor was 2.5, indicating a substantial amount of excess 
scatter beyond the internal error estimates. We do not 
believe this scatter to be due to the effect of starspots, 
given that the Hipparcos photometry of this system in- 
dicates a scatter of no more than 5 mmag; the resulting 
maximum starsp ot-induced astrometric noi se would be 
~ 4/iarcseconds (jMuterspaugh et aLll2006al ). We have 
however identified possible instrumental sources of this 
sys tematic error and develope d methods for reducing it, 
see IMuterspaugh et aTl (|2007[ ) . Nevertheless, the exist- 
ing astrometry is sufficient to detect astrometric motion 
induced by the short-period subsystems. 

2.2. Previous Astrometry 

In addition to our astrometry, 88 Tau A has 
been followed by a number of observers with speckle- 
interferometric techniques. We use 20 observations tabu- 
lated in the 4th Catalog of Interferometric Measurements 
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Fig. 2.— Two spectra of 88 Tau, (a) JD 2,452,539.97 and (b) JD 
2,454,006.01, compared with (c) the spectrum of the lAU radial 
velocity standard 10 Tau. The four components of 88 Tau A are 
identified for two different Unes. Solid tick marks indicate Unes of 
the Aa (3.57 day period) binary with the longer tick mark iden- 
tifying the primary. Dotted tick marks indicate lines of the Ab 
(7.89 day period) binary with the longer tick mark identifying the 
primary. 

of Binary Stars^ (|Hartkopf et al.ir200lD to further con- 
strain our fit. Although of somewhat lower precision, the 
considerable time-baseline (including observations dat- 
ing from 1985) help constrain the parameters of the wide 
orbit. In many cases the published astrometry lacks un- 
certainties, and we therefore assigned a plausible initial 
uncertainty of 3 mas in separation and 2 degrees in posi- 
tion angle to these points (we used published uncertain- 
ties for the points where such were available). We then 
performed a least-squares fit of a single Keplerian orbital 
model (corresponding to the Aa-Ab orbit - the subsys- 
tems are far too small to be detected by these data), and 
scaled all of the uncertainties so as to yield a reduced 
of unity. We find the average uncertainty in separation 
to be 5 milli-arcseconds, and the average position-angle 
uncertainty to be 3 degrees. 

2.3. Spectroscopic observations and reductions 

From 1984 January through 2006 September we ob- 
tained 82 spectrograms of 88 Tau with the Kitt Peak Na- 
tional Observatory (KPNO) 0.9 m coude feed telescope, 
coude spectrograph, and a TI CCD detector. Sixty-eight 
spectrograms are centered in the red at 6430 A, cover a 
wavelength range of about 80 A, and have a two pixel 
resolution of 0.21 A. Those spectra have signal-to-noise 
ratios of ~250. The remaining 14 spectrograms are cen- 
tered in the blue at 4500 A, cover a wavelength range 
of 85 A, and have a resolution of 0.22 A. Signal-to-noise 
ratios of ^^300 are typical. 

From 2004 January through 2007 April we acquired 
29 spectrograms with the Tennessee State University 
2 m automatic spectroscopic telescope (AST), fiber- fed 
echelle spectrograph, and a 2048 x 4096 SITe ST-002A 
CCD. The echelle spectrograms have 21 orders, cover- 
ing the wavelength range 4920-7100 A with an average 
resol ution of 0.17 A. The typical signal-to-noise ratio is 
~50. lEaton k WilliamsonI (j20M) have given a more ex- 
tensive description of the telescope, situated at Fairborn 
Observatory near Washington Camp in the Patagonia 

* |http : //ad . usno ■ navy . mll/wds/lnt4 . html | 
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TABLE 1 

PHASES ASTROMETRIC DATA FOR 88 TAU A 



[JD-2400000.5 


ARAcos(5) 


ADec 






4^e 




0"I3ec 




N 


ADC 


Align 


Rate 


rmn 


maj 








(iiia,s) 












t 










52979.34108 


-32.1751 


-101.5913 


20.6 


477.9 


163.71 


458.8 


135.5 


-0.98733 


2353 








100 


53034.13365 


-40.3179 


-91.5425 


14.0 


311.9 


152.78 


277.4 


143.2 


-0.99394 


2951 








100 


53250.50169 


-65.6825 


-53.5438 


40.2 


1775.4 


147.24 


1493.2 


961.3 


-0.99876 


1027 








100 


53271.46723 


-70.0370 


-48.2688 


16.6 


669.7 


150.15 


580.9 


333.7 


-0.99836 


4005 








100 


53291.40191 


-71.5816 


-44.8765 


88.3 


1759.4 


148.37 


1498.7 


925.8 


-0.99371 


408 








100 


53294.47661 


-71.0249 


-44.6002 


32.3 


1700.7 


164.09 


1635.6 


467.2 


-0.99742 


1355 








100 


53312.38113 


-74.1285 


-40.9475 


7.8 


68.0 


154.15 


61.3 


30.5 


-0.95879 


8150 








100 


53320.33408 


-74.5082 


-39.6052 


40.3 


2076.1 


150.23 


1802.1 


1031.5 


-0.99899 


1838 








100 


53340.29304 


-77.4848 


-35.5088 


18.2 


215.1 


152.98 


191.8 


99.1 


-0.97862 


3602 








100 


53341.28228 


-77.5104 


-35.2746 


17.2 


551.5 


150.70 


481.0 


270.3 


-0.99734 


3630 








100 


53605.52764 


-104.7426 


15.6166 


18.3 


956.5 


147.64 


808.1 


512.1 


-0.99911 


2855 








100 


53606.51559 


-106.5080 


16.8975 


15.7 


972.4 


146.14 


807.5 


542.0 


-0.99939 


3205 








100 


53614.50698 


-105.6799 


17.1892 


25.4 


773.7 


147.94 


655.8 


411.2 


-0.99735 


1716 








100 


53687.37864 


-112.2999 


31.2082 


24.7 


227.4 


31.80 


193.7 


121.7 


0.97114 


3260 








100 


53711.29715 


-114.6779 


35.8970 


80.7 


880.5 


29.41 


768.1 


438.0 


0.97743 


833 


1 





100 


53712.28437 


-114.0278 


36.0764 


41.2 


326.1 


27.87 


288.9 


156.8 


0.95491 


2889 


1 





100 


53789.14419 


-120.0382 


51.3801 


59.4 


3107.9 


40.17 


2375.3 


2005.2 


0.99925 


761 


1 





100 


53790.13775 


-119.0373 


52.0709 


77.5 


2663.3 


39.39 


2058.8 


1691.3 


0.99824 


624 


1 





100 


54030.46119 


-135.6408 


94.2582 


50.4 


1514.6 


163.75 


1454.1 


426.6 


-0.99240 


1268 


1 


1 


50 


54055.38631 


-136.2430 


98.5452 


48.0 


870.1 


161.74 


826.4 


276.5 


-0.98318 


1679 


1 


1 


50 


54061.37247 


-136.6968 


99.4972 


16.3 


302.1 


163.08 


289.0 


89.3 


-0.98165 


4216 


1 


1 


50 


54075.33188 


-136.8574 


101.6204 


13.6 


250.1 


162.18 


238.2 


77.7 


-0.98300 


5902 


1 


1 


50 


54083.30564 


-137.8524 


103.1506 


31.6 


688.5 


161.39 


652.6 


221.8 


-0.98863 


3845 


1 


1 


50 


54084.31595 


-136.7985 


103.1084 


63.9 


1741.5 


163.11 


1666.5 


509.6 


-0.99139 


1182 


1 


1 


50 


54103.25456 


-137.6742 


106.3338 


42.7 


1011.9 


162.34 


964.3 


309.6 


-0.98949 


743 


1 


1 


50 


54138.15348 


-139.3256 


112.2494 


28.6 


457.9 


161.38 


434.1 


148.7 


-0.97916 


3142 


1 


1 


50 



Note. — All quantities are in the ICRS 2000.0 reference frame. The uncertainty values presented in these data have been scaled 
by a factor of 2.5 over the formal (internal) uncertainties for each night. Column 6 {4'e) is the angle between the major axis of the 
uncertainty ellipse and the right ascension axis, measured from increasing differential right ascension through increasing differential 
declination. N is the number of scans obtained in a night; each scan typically represents 0.5-1 second of integration. AE)C indicates 
that the observations made use of the automatic dispersion compensator. Rate indicates the tracking rate of the fringe tracker 
used to stabilize the fringe phase during measurement. Align indicates whether or not the automatic alignment system was used 
to stabilize the system pupil. 



Mountains of southeastern Arizona, and its operation. 

For the KPNO spectrograms we determined radial ve- 
locities wn^h_yieIRAF cross-correlation program FX- 
COR (jFitzpatrickl Il993f) . fitting Gaussian functions to 
the individual cross-correlation peaks. Double Gaussian 
fits were required to obtain individual velocities from 
blended cross-correlation peaks. The lAU radial velocity 
standard star 10 Tau was used as the cross-correlation 
reference star for the red-wavelength spect rograms. Its 
veloci ty of 27.9 km s~^ was adopted from IScarfe et all 
()1990() . Lines in the wavelength region redward of 6445 A 
are not particularly suitable for measurement because 
most features are blends, and there are a number of 
modest strength water vapor lines. Thus, the radial ve- 
locities were determined from lines in the region 6385- 
6445 A. However, this 60 A portion of the spectrum is so 
small that a spectrum mismatch, caused by the varying 
strength of line blends with temperature, between the A 
an d F spectral type compon ents of 88 Tau and the F9 IV- 
vdl eenan fc McNeillll989f ) spectral type of the reference 
star 10 Tau, can significantly alter the measured velocity. 
Thus, instead of cross-correlating this entire 60 A wave- 
length region, only the wavelength regions around two 
or three of the strongest and least-blended lines, usually 
the Fe I lines at 6394 and 6412 A plus the Ca I line at 
6439 A, were cross-correlated. 

At blue wavelengths the Am star dominates the spec- 
trum. To compute velociti es from those spectro grams, 
68 Tau, spectral type A2 IV (jAbt fc Morrellll 1 995( 1 . which 



has a velocity of 39.0 km s'^ (|Fekellll999l ). was used as 
the reference star. The region between 4485 and 4525 A 
was cross correlated. Velocities for our KPNO spectra 
are given in Table [2] 

For the Fairborn Observatory AST spectra, lines in 
approxir nately 100 regions, centered on the rest wave- 
lengths (|Moore et al.l |1966[) of relatively strong lines 
(mostly of Fe I and Fe II) that were not strong blends, 
were measured. Lines at the ends of each echelle order 
were excluded because of their lower signal-to-noise ra- 
tios. A Gaussian function was fitted to the profile of 
each component. Double Gaussian fits were required to 
represent blended components. The difference between 
the ob served wavelength a nd that given in the solar line 
list of lMoore et al.l ()1966( l was used to compute the ra- 
dial velocity, and a heliocentric correction was applied. 
The final mean velocity for each observation is given in 
Table [H Unpublished velocities of several lAU standard 
stars with F dwarf spectral types indicate that the Fair- 
born Observatory velocities have a small zer o-point offset 
of —0 .3 km s~^ relative to the velocities of IScarfe et al.l 
(fT99(ih . 

2.4. Preliminary spectroscopic analysis 

In a study of lithium in Am stars iBurkhart fc Couprvl 
(198^ acquired two high-resolution spectrograms of 88 
Tau in the 6710 A region. Comparing two sets of lines 
in the two spectra, they reported detecting the lines of 
5 different components. Figure 2 presents two spec- 
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TABLE 2 

KPNO Radial Velocity data for 88 Tau A 



HJD-2400000.5 


VAal 


WeigW 


VAa2 


Weight'' 


VAbl 


Weighf^ 


VAb2 


Weight<i 




(kms-1) 




(kms-1) 




(kms-1) 




(kms-1) 




45718.327 


91.4 


1.0 


-66.5 


1.0 


- 


- 


- 


- 


46388.389 


-48.7 


1.0 


141.6 


1.0 


40.1 


1.0 


1.1 


1.0 


46390.212 


99.8 


1.0 


-85.2 


1.0 


-5.8 


1.0 


- 


- 


46718.473 


96.5 


1.0 


-93.7 


1.0 


59.8 


0.4 


-4.5 


1.0 


46720.487 


-51.6 


1.0 


139.5 


1.0 


- 


- 


- 


- 


47152.190 


-45.3 


1.0 


123.6 


1.0 


57.9 


1.0 


- 


- 


47245.112 


-50.2 


1.0 


131.5 


1.0 


- 


- 


- 


- 


47456.264 


-52.2 


1.0 


129.7 


1.0 


9.3 


1.0 


- 


- 


47556.098 


-57.7 


1.0 


135.8 


1.0 


- 


- 


- 


- 


47624.100 


-53.5 


1.0 


125.7 


1.0 


- 


- 


- 


- 


47626.145 


66.5 


1.0 


-57.7 


1.0 


- 


- 


4.3 


1.0 


47627.123 


-49.3 


1.0 


118.9 


1.0 


- 


- 


- 


- 


48345.123 


-60.5 


1.0 


133.2 


1.0 


- 


- 


- 


- 


48347.103 


96.5 


1.0 


-102.4 


1.0 


9.4 


1.0 


69.9 


0.4 


48356.100 


-58.6 


1.0 


128.6 


1.0 


10.0 


1.0 


59.1 


1.0 


48505.518 


-46.5 


1.0 


106.1 


1.0 


14.0 


1.0 


53.0 


1.0 


48573.407 


-49.1 


1.0 


112.7 


1.0 


42.5 


0.4 


26.2 


0.4 


48604.284 


90.2 


1.0 


-98.1 


1.0 


57.7 


0.4 


15.4 


1.0 


48607.306 


75.8 


1.0 


-74.8 


1.0 


3.0 


1.0 


- 


- 


48913.412 


-43.0 


1.0 


105.6 


1.0 


19.0 


1.0 


50.9 


1.0 


48916.375 


-55.5 


1.0 


125.2 


1.0 


- 


- 


- 


- 


49246.472 


76.6 


1.0 


-67.1 


1.0 


6.6 


0.4 


- 


- 


49248.525 


-56.1 


1.0 


131.7 


1.0 


44.0 


1.0 


15.3 


1.0 


49250.537 


101.6 


1.0 


-93.6 


1.0 


58.9 


0.4 


1.7 


1.0 


49302.418 


-53.7 


1.0 


132.4 


1.0 


12.9 


1.0 


51.9 


1.0 


49307.374 


91.3 


1.0 


-87.2 


1.0 


23.0 


0.4 


34.6 


0.4 


49618.436 


94.0 


1.0 


-92.1 


1.0 


- 


- 


- 


- 


49622.425 


60.6 


1.0 


-41.1 


1.0 


- 


- 


- 


- 


49677.381 


-50.0 


1.0 


132.6 


1.0 


- 


- 


- 


- 


49971.535 


93.8 


1.0 


-78.7 


1.0 


-9.8 


1.0 


53.5 


0.4 


49973.434 


-48.2 


1.0 


138.8 


1.0 


- 


- 


- 


- 


49973.535 


-51.1 


1.0 


142.9 


1.0 


17.0 


0.4 


26.0 


0.4 


50364.509 


102.5 


1.0 


-85.3 


1.0 


5.8 


0.4 


30.0 


0.4 


50366.371 


-48.5 


1.0 


146.0 


1.0 


- 


- 


- 


- 


50400.392 


100.7 


1.0 


-88.4 


1.0 


38.2 


1.0 


-1.3 


1.0 


50404.270 


91.1 


1.0 


-66.5 


1.0 


- 


- 


- 


- 


50721.423 


93.8 


1.0 


-67.1 


1.0 


-9.9 


1.0 


38.0 


0.4 


50721.501 


96.1 


1.0 


-77.1 


1.0 


-11.3 


1.0 


40.2 


1.0 


50755.460 


-46.1 


1.0 


139.7 


1.0 


38.8 


1.0 


-11.3 


0.4 


50757.384 


104.0 


1.0 


-86.7 


1.0 


34.0 


1.0 


-2.7 


1.0 


50832.293 


104.7 


1.0 


-82.2 


1.0 


-9.7 


1.0 


44.2 


1.0 


50833.156 


57.7 


1.0 


- 


- 


- 


- 


- 


- 


51088.384 


-7.6 


1.0 


92.2 


1.0 


- 


- 


- 


- 


51089.481 


106.0 


1.0 


-85.2 


1.0 


- 


- 


- 


- 


51091.427 


-49.1 


1.0 


149.9 


1.0 


-18.8 


0.4 


47.3 


1.0 


51093.343 


103.8 


1.0 


-79.4 


1.0 


3.6 


0.4 


22.0 


0.4 


51473.377 


-51.8 


1.0 


147.2 


1.0 


43.5 


1.0 


-18.3 


0.4 


51475.281 


107.5 


1.0 


-85.7 


1.0 


- 


- 


- 


- 


51475.362 


107.8 


1.0 


-85.7 


1.0 


22.1 


0.4 


4.1 


0.4 


51803.489 


95.8 


1.0 


-66.7 


1.0 


- 


- 


- 


- 


51805.463 


-49.1 


1.0 


147.2 


1.0 


41.6 


1.0 


-21.2 


0.4 


51807.463 


107.0 


1.0 


-86.9 


1.0 


5.1 


0.4 


22.7 


0.4 


52016.110 


-46.8 


1.0 


142.6 


1.0 


- 


- 


- 


- 


52180.462 


-49.0 


1.0 


146.2 


1.0 


-14.7 


1.0 


43.8 


1.0 


52182.481 


104.0 


1.0 


-86.2 


1.0 


28.2 


0.4 


0.4 


0.4 


52327.104 


-48.8 


1.0 


145.5 


1.0 


35.9 


1.0 


-0.9 


0.4 


52329.148 


89.8 


1.0 


-64.0 


1.0 


-8.0 


0.4 


38.2 


0.4 


52537.494 


-50.0 


1.0 


137.7 


1.0 


35.0 


1.0 


-4.0 


0.4 


52539.468 


106.0 


1.0 


-89.1 


1.0 


42.9 


1.0 


-7.0 


1.0 


52541.378 


-46.3 


1.0 


144.8 


1.0 


- 


- 


- 


- 


52541.473 


-45.5 


1.0 


139.3 


1.0 


1.5 


0.4 


38.9 


1.0 


52705.198 


-45.8 


1.0 


126.6 


1.0 


48.5 


1.0 






52707.190 


101.6 


1.0 


-88.8 


1.0 










52709.127 


-54.6 


1.0 


144.6 


1.0 






51.1 


1.0 


52903.438 


93.6 


1.0 


-75.7 


1.0 










52941.345 


-50.8 


1.0 


141.1 


1.0 


51.1 


1.0 


-6.8 


0.4 


53273.479 


-53.1 


1.0 


139.3 


1.0 


41.0 


1.0 


7.9 


1.0 


53278.490 


93.6 


1.0 


-84.9 


1.0 










53637.454 


-55.3 


1.0 


134.5 


1.0 


12.9 


1.0 


42.4 


1.0 


54001.454 


-42.5 


1.0 


113.7 


1.0 


3.6 


0.4 


66.6 


1.0 


54003.423 


89.9 


1.0 


-90.3 


1.0 










54005.351 


-58.5 


1.0 


137.5 


1.0 


62.1 


1.0 






54005.506 


-57.7 


1.0 


134.6 


1.0 


60.4 


1.0 


0.9 


0.4 



^ An observation of Aal of unit weight has a standard error of 2.0 kms ^.'^ An observation of Aa2 of unit weight 
has a standard error of 2.3 kms"^.*^ An observation of Abl of unit weight has a standard error of 2.6 kms~^.*^ An 
observation of Ab2 of unit weight has a standard error of 2.6 kms~^. 
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TABLE 3 

Fairborn Observatory Radial Velocity data for 88 Tau A 



HJD-2400000.5 


VAal 


WeigW 


VAa2 


Weight'' 




Weighf^ 


VAb2 


Weigh 




(kms-1) 




(kms-1) 




(kms-1) 




(kms-1) 




53020.204 


-38.5 


1.0 


118.0 


1.0 


50.8 


1.0 






53032.197 


101.1 


1.0 


-90.3 


1.0 


-6.3 


1.0 


53.7 


1.0 


53052.273 


-43.0 


1.0 


126.0 


1.0 


50.9 


1.0 






53276.489 


-37.3 


1.0 


110.0 


1.0 






53.0 


1.0 


53285.468 


85.9 


1.0 


-69.9 


1.0 


10.8 


1.0 


41.8 


1.0 


53314.392 


102.0 


1.0 


-95.0 


1.0 


8.8 


1.0 


40.6 


1.0 


53350.390 


93.6 


1.0 


-85.3 


1.0 


46.3 


1.0 


7.7 


1.0 


53395.275 


-28.2 


1.0 


97.3 


1.0 


-2.4 


1.0 


58.0 


1.0 


53405.184 


-44.4 


1.0 


127.0 


1.0 


39.9 


1.0 


11.9 


1.0 


5.36.30.518 


-55.6 


1.0 


139.5 


1.0 


2.7 


1.0 


57.9 


1.0 


53644.477 


-48.0 


1.0 


127.2 


1.0 


38.3 


1.0 


20.1 


1.0 


53659.486 


-33.5 


1.0 


102.9 


1.0 


58.4 


1.0 






53700.387 


90.5 


1.0 


-82.3 


1.0 










53741.243 


-56.5 


1.0 


139.1 


1.0 


-1.2 


1.0 


56.7 


1.0 


54191.132 


-61.0 


1.0 


136.8 


1.0 


7.3 


1.0 


63.7 


1.0 


54194.134 


-28.7 


1.0 


85.1 


1.0 










54194.156 


-30.2 


1.0 


89.5 


1.0 










54198.110 


-56.9 


1.0 


132.4 


1.0 











^ An observation of Aal of unit weight has a standard error of 2.0 kms ^.'^ An observation of Aa2 of unit weight 
has a standard error of 2.3 kms"^.*^ An observation of Abl of unit weight has a standard error of 2.6 kms-^.*^ An 
observation of Ab2 of unit weight has a standard error of 2.6 kms~^. 
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tra of 88 Tau A in the 6430 A region that show the 
two components of the 3.57 day binary near opposite 
nodes in their spectroscopic orbit, when the components 
have their maximum velocity separation. Between the 
two "outside" hues are two additional weak components. 
From a careful inspection of our KPNO spectra, as well 
as the ones obtained at Fairborn Observatory, we find 
lines of only 4 components rather than the 5 reported by 
iBurkhart fc Couprvl (|1988[ ). In Figure [H many lines of 
the Am star are ^--^5% deep, while the lines of the other 3 
components typically have line depths <2.5%. Correctly 
identifying components in such a weak-lined and com- 
plex spectrum is not easy because in many spectra two 
or more of the components are blended. 

While the identification of the components of the 3.57 
day binary is straightforward, the very weak lines of the 
other two components are similar in strength and line 
width, making it difficult to tell them apart. To deter- 
mine a preliminary orbital period, we initially examined 
only the latter portion of our KPNO velocities, obtained 
from MJD 50400 to 54000. For each observation we com- 
puted the absolute value of the velocity difference be- 
tween the two components and then used those results 
as the input data for two different period finding ap- 
proaches. First, a sine curve was fitted to the velocity 
differences for trial periods between 1 and 100 days with 
a step size of 0.0005 days. The period with the small- 
est sum of the squared residuals was adopted as the best 
period. Next, a search over a similar period range and 
with the same step size was done with the least string 
method (Bopp et al. 1970). Both searches resulted in a 
period of 3.9435 days. Doubling this period produced an 
orbital period of 7.887 days. Separate analyses of our 
earlier KPNO velocities as well as the Fairborn Observa- 
tory velocities produced a similar orbital period. We then 
adopted the 7.887 day period and computed a phase dia- 
gram to identify correctly the components. Afterward we 
compared those results with an attempt at visual iden- 
tification, based on which set of lines appeared to be 
stronger in each spectrum. The visual inspection cor- 
rectly identified the more massive component only about 
half of the time. Apparently, the lines of these two com- 
ponents are similar enough that weak lines from other 
components and noise can significantly affect the appar- 
ent line strengths. Thus, we conclude that in our spec- 
tra it is not possible to differentiate the two components 
based on line strength. 

2.5. Orbital Models 

In modeling the hierarchical quadruple system we 
make the simplifying assumption that the three orbital 
systems do not perturb each other during the time of 
our observations, i.e. we use three Keplerian orbital sys- 
tems, one wide (Aa-Ab) and slow (18-year period), and 
two short period systems: Aal-Aa2, 3.57-day period, and 
Abl-Ab2, 7.89-day period. Note that one cannot simply 
superimpose the separation vectors from the three mod- 
els; this is because the PHASES observable is the angle 
between the two Centers-of-Light (COL) of the short- 
period systems. 

J/obs = Aa-Ab 

+ (l+i?Aj(l+i^,)'^^^-^^' 



-{l + RAb){l + LAbf^''-^'" 
Here RAa = -^^Aa2/-^Aai IS the Aa component mass ra- 
tio and Lau — -^Aa2/iAai the luminosity ratio, while 
RAb = MAb2/MAbi and LAb = -^Ab2/iAbi are the corre- 
sponding ratios for the Ba-Bb sub-system. Including this 
coupling term for astrometric data is important when a 
full analysis, including radial velocity data, is made. 

3. RESULTS 

The best-fit orbital model was found with an iterative 
non-linear least-squares minimization scheme. The best- 
fit parameters are found in Table ID The reduced Xr of 
the combined fit to PHASES, radial velocity, and previ- 
ous differential astrometry data is 1.37. This combined 
set has 378 data points (49 of which are two-dimensional 
astrometric points) and 23 free parameters. The fits to 
the astrometric and radial-velocity data for the various 
subsystems are shown in Figures 3-7. We find that the 
two short-period systems have eccentricities consistent 
with zero, and we therefore held these parameters fixed 
at zero for the fit. The time of maximum primary appar- 
ent velocity is chosen as zero orbital phase. To investi- 
gate the consistency of two astrometric data sets we also 
re-ran the fit, without including the previous astrome- 
try. We find that the final results are consistent to within 
one sigma; including the previous astrometry does reduce 
certain parameter uncertainties by a small amount. In 
Table [4] we compare ou r best-fit values with those avail- 
able in the literature (|Abt fc Levvl 119851 : iBalega et al.l 
I1999D where possible; we find generally good agreement. 

Finally, we considered the possibility of an additional 
massive body in the system. At this point, a full mass- 
period phase space search for a fifth component is com- 
putationally prohibitive and beyond the scope of this in- 
vestigation. However, through the course of the inves- 
tigation, we discovered a possible periodicity that war- 
ranted follow-up. Before we had determined the 7.89 
day period of the Abl-Ab2 system, we had searched the 
PHASES data for this third orbit signal with the period- 
searching program used to find the signature s of addi- 
tional companions iMuterspaugh et al.l l2006c[ ): because 
the Aa-Ab and Aal-Aa2 periods were known already, 
these orbits were seeded at their best fit double Kepler 
values. While radial velocities eventually isolated the pe- 
riod for Abl-Ab2 as 7.89 days, this astrometric search 
additionally identified a curious improvement near 56 
days. The size of that orbit could have corresponded 
either to a brown dwarf or to a nearly equal luminosity 
star-star subsystem. Inspired by this identification, a 5- 
component, quadruple Keplerian model was seeded with 
the three known orbits (Aa-Ab, Aal-Aa2, and Abl-Ab2) 
and a fourth perturbation with a period of 56 days. Sep- 
arately, the harmonic of 28 days was also investigated. 
This model included all astrometric and radial velocity 
data. The 5-component model with a candidate period 
near 28 days (C/Ab mass ratio 0.044) yields a signif- 
icant improvement in the reduced of the fit over a 
4-body model (491 for 30 parameters and 397 degrees of 
freedom vs. 547 for 23 parameters and 404 degrees of 
freedom, corresponding to a ~ 2.3a improvement from 
the ratio test). The candidate period near 56 days 
yields a best-fit ~ 474 for a 3.7cr improvement over 
the 4-component model, yet with a very large best-fit 
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TABLE 4 
Best-fit Orbital Parameters for 



Tau a 



TABLE 5 

Derived System Parameters for. 88 Tau A 



Parameter 


Value & Uncertainty 


Previous Value 


, 2 
X 


547 




X.r 

No. Par am. 


1.37 




zo 




PAaAb(<iays) 




uoyo m 4^ 


^AaAb 


f) nvi -1- n nn2fi 


n (183 + OOS'^ 

KJ.yJiJ'-J _1_ W.L/VJCJ 


1 * n -1 ( r\ oo" 1 
^AaAbV*-'^'-&-7 


KO qo'^ + n (148 


70.4 lb 0.4^^ 


(■(J A A 1_ I Hp0' 1 

"^AaAbV'-^'-&'/ 


205 7+ 1 2 


222.0 lb 3.0*^ 


AaAb \LLL\i.oi^ ) 


55261 zb 22 




O A ^ A K fdee;.) 

AaA D V 6 / 


146.734 ± 0.067 


146.6 lb 0.3*= 


MAa(MQ) 


3.42 ± 0.18 




MAhfMm) 


2.13 ± 0.13 






50.70 lb 0.88 




PAalAa2 (days) 


Q t;7i oQfi -1- n nnnnn'i 
o.oi -i-uyu ^ u.uuuuuo 


Q t^yi QQ1 _|_ fi nnnnoiR^ 
o.o/ioyj- ^ u.uuuuuo 


eAalAa2 


U.U 


U 


iAalAa2 (deg.) 






"^AalAa2 V-^^B v 


0.0 




^ 0,AalAa2 l^iVinju; 




OS 1 (I'i -1- (1 01 


''AalAaz } 






MAa2/MAal 


0.6602 ± 0.0028 




LAa2/LAai (K-band)=' 


0.249 ±0.035 




PAbiAb2(days) 


7.886969 ± 0.000066 




eAblAb2 


0.0 




iAblAb2(deg.) 


27.23 ±0.72 




'^AblAb2(deg.) 


0.0 




To,AbiAb2{MHJD) 


52507.31 ±0.02 




r2AbiAb2(deg.)'' 


34.0 ±8.2 




MAb2/MAbl 


0.988 ± 0.024 




LAb2/LAbi (K-band)'' 


0.865 ±0.028 




Vo{KPNO,kms-i) 


23.70 ±0.17 




Vo{Fairborn, kms~^) 


23.91 ±0.31 





Parameter 



eccentricity (0.93). However, we note that the 28-day 
period is close to the product of the Aa and Ab orbital 
periods (~ 28.1 days) and the 56 day period is a har- 
monic of this, raising suspicion regarding the reliability 
of the solution. Hence the astrophysical significance of 
this additional component remains unclear. We do not 
claim this as a detection, but it is an intriguing result 
worthy of observational follow-up. 

3.1. Relative Orbital Inclinations 
The mutual inclination $ of two orbits is given by 

cos $ — cos ii cos 12 + sin ii sin 12 cos (fii — £^2) (2) 

where ii and 12 are the orbital inclinations and fli 
and ri2 are the longitudes of the ascending nodes. For 
this quadruple system we derive three separate mutual 
inclinations, corresponding to the three possible pair- 
wise comparisons of the three orbits in this system (i.e. 
Aab - AalAa2, Aab ~ AblAb2 and AalAa2 - AblAb2). 
The resulting values found from our combined orbital 
solutions are given in Table [S] The large mutual inclina- 
tions indicate that the system orbits are not co-planar. 

It should be noted that even with both COL- 
astrometry and radial-velocity data, there exists a pa- 
rameter degeneracy corresponding to an exchange of 
the ascending and descending nodes together with a 
change in the luminosity ratio (interchanging which is 
the brighter star). Given one solution for the mass and 
luminosty ratios {R and Li), the other possible luminos- 



*AaAb-AalAa2 (deg.) 
*AaAb-AblAb2(deg.) 

*AalAa2-AblAb2 (deg.) 
7r(asec) 

aAaAb (milliarcsec) 
aAalAa2 (miUiarcsec) 
aAalAa2,col (miUiarcsec) 
aAblAb2 (milliarcsec) 
aAblAb2 col (miUiarcsec) 
aAB(AU) 

aAalAa2 (AU) 

aAblAb2(AU) 

MAal(M0) 

MAa2(M0) 

MAbl(M0) 

MAb2(Mo) 

MK,Aal 

MK,Aa2 

MK,Abl 



Value 


Uncertainty 


143.3 


±2.5 


82.0'' 


±3.3 


115.8*= 


±4.6 


0.01973 


±0.00034 


240.1 


±5.3 


1.359 


±0.034 


0.270 


±0.032 


1.967 


±0.054 


0.065 


±0.020 


12.17 


±0.17 


0.0689 


±0.0012 


0.0997 


±0.0021 


2.06 


±0.11 


1.361 


±0.073 


1.069 


±0.069 


1.057 


±0.068 


0.69 


±0.26 


2.20 


±0.28 


2.31 


±0.27 


2.00 


±0.27 



Note. — The parameters derived from the best- 
fit model values in Table [4] and their uncertainties. 
"COL" refers to the amplitude of the motion of the 
Center of Light of the subsystem in question. Note 
that our combined astrometry and radial velocity 
model fits for the system masses directly. 
^ An alternate solution has ^AaAb — AbiAb2 — 58 
deg. if LAb2/LAbi > 1-^ An alternate solution has 

^AalAa2-AblAb2 ^ 107 dcg. if L^taZ-^-Afal > 1- 



An alternate, but disfavored, solution has LAa2/LAai — 1.48 and 
AaiAa2 — 326.^ Au alternate solution has LAb2/LAbi — 1.10 and 
AbiAb2 ^ 205.'' From fBalega et al.l (fT999h '^ From [Abr&: Levvl ||1985D 



ity ratio can be found from 

2R + RLi 



\^2Li-R 



(3) 



In a quadruple stellar system such as 88 Tau A there are 4 
possible model solutions. However, as can be seen in Fig. 
[2] the luminosity of the Aa\ component is clearly greater 
than the Aa2 component, hence we choose the solution 
where LAaiAa2 = 0.24. However, given the nearly-equal 
masses of the Abl and Ab2 components, it is not entirely 
clear which is the preferred solution {LAbiAb2 = 0.87 or 
1.10), and we calculate the two possible values for the 
corresponding mutual inclinations. 

The two possible mutual inclination values of the Ab 
system are both in the range for inc lination-ecce ntricity 
oscillations of 39.2 - 140.8 degrees ("Kozail 11963 ). while 
the Aa system is not (though only 1 standard deviation 
away from the limit). During these "Kozai-cycles" the or- 
bital eccentricity varies on a timescale oc Pg^^/Pmi^ 10^ 
years for the Ab system; Kiseleva-Eggleton & Eggleton 
2001), and in the absence of damping factors the eccen- 
tricity values would range from ~ to ^ 0.98 (Kiseleva 
et al., 1998; however, see Kiseleva-Eggleton & Eggle- 
ton 2001 for a discussion of various factors that may 
limit these excursions). A full treatment of the dynam- 
ics of this system, including the effects of tidal friction, 
quadrupolar distortion and general relativity is beyond 
the scope o f this paper. 

Recently, iFabrvckv fc Tremaiii3 (|2007f l studied the ef- 
fect that Kozai oscillations would have on the distribu- 
tions of orbital properties of triple systems. Systems 
with mutual inclinations in the range where Kozai cy- 
cles occur evolve rapidly as tidal dissipation during the 
high-eccentricity (and hence close approach) phase of the 
oscillation causes the orbit of the inner binary to shrink. 
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88 Tau A Visual Orbit 
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Fig. 3. — (left)The best-fit visual orbit of the 88 Tau Aa-Ab system, together with previously available astrometric data, and our PHASES 
astrometry. Note that the error ellipses of the PHASES data appear smaller than the points used to indicate the data, (right) A close-in 
view of a subsection of the PHASES astrometry, together with the best-fit orbital model. 



Aa1-Aa2 Orbit (145 deg., Aa-Ab, Ab1-Ab2 Subtracted) 



Aa1 -Aa2 Orbit (38 deg., Aa-Ab, Abl -Ab2 Subtracted) 
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Fig. 4. — The astrometric orbit of the 88 Tau Aal-Aa2 subsystem, projected along two different axes (145 degrees East of North deg on 
the left, 38 degrees on the right.) In each case the motion of the Aa-Ab and other sub-system has been removed. The axes correspond to 
the two most common orientations of the minor axis of the positional error ellipses (which vary slightly from night to night, and between 
baselines). For clarity, only those observations where the projected uncertainty is less than 300 fias have been includeci in the plot (all 
observations are included in the fit.) 

The resulting mutual inclination distribution is strongly 
peaked near the critical values of 39 and 141 degrees. It 
is therefore interesting to note that the Aa system is so 
close to one of these predicted values. The Ab system is 
more ambiguous: the preferred (albeit only weakly) so- 
lution for the Ab system yields a mutual inclination that 
would be expected to result in rapid orbital evolution. In- 
terestingly, the second possible solution is also relatively 
close to the critical limit. Clearly, further observations 
will be needed in order to remove the ambiguity in the 
node and hence mutual inclination. 

3.2. Component Masses and Distance 

This study represents the first determination of the 
orbital inclinations in this system, and hence the first 
time the masses have been determined; the precision 



achieved is ~ 5% for the Aa components and ~ 6 % 
for the Ab components. The parallax is found to be 
19.73 ± 0.34 mas (1.7% uncertainty), placing the system 
at a greater distance than that estimated by Hipparcos 
(21.68 ± 0.82 mas; ESA 1997). Our greater distance re- 
solves the nia,ss/lu minosity discrepancy pointed out by 
iBalega et all (|1999[ ). which arises if one assumes the Hip- 

p arcos distance to this s ystem. 

iSchoeller et"an (|1998[ ) and IBalega et all (|1999[ ) stated 
that the 3.57 day binary is eclipsing but gave no ref- 
erence for this claim. Our inclination for that binary 
is 110.6°, a value that is about 10° above the maximum 
value required for eclipses to occur if c anonical va lues are 
assumed for the radii of the two stars (|Gravlll99^ . Thus, 
the system is not eclipsing. We also note that because of 
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Ab1-Ab2 Orbit (145 deg., Aa-Ab, Aa1-Aa2 Subtracted) 
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Fig. 5. — The astrometric orbit of the 88 Tau Abl-Ab2 sub- 
system. The motion of the Aa-Ab and other sub-system has been 
removed. The separations shown are projected along an axis ori- 
ented 145 degrees East of North, corresponding to mean orientation 
of the minor axis of the positional error ellipses (which vary slight 
from night to night). For clarity, only those observations where the 
projected uncertainty is less than 300 /las have been included in 
the plot (all observations are included in the fit). 

our lower inclination, o ur masses are about 0.2Mq larger 
than those adopted by iBalega et al.l ([IQQQT ). 

3.3. Component Luminosities 

As part of the combined astrometric and radial veloc- 
ity fit we can solve for the K-band luminosity ratios of 
the components; this is because the distance and sub- 
system total masses are essentially determined by the 
observations of the wide Aa-Ab system, while the sub- 
system mass ratios are found from the subsystem radial 
velocities; this leaves only the component luminosity ra- 
tios dependent on the size of the observed astrometric 
perturbation. 

While PTI cannot provide precise determinations of 
the total system magnitude mx or the Aa-Ab system 
differentia l magnitude Awj<-, these can be found in the 
literature. iBalega et all (|2001h gave Atuk' = 1.29 ±0.12 
for the Aa-Ab system. We derive the resulting absolute 
K magnitudes and list them in Table [SJ the results are 
dominated by the uncertainty in the total magnitude. 
Note that we assume the solution where LAb2/LAbi < 1- 
We compare the determined masses and K-band ab- 
solute magnitudes to published theoretical isochrones 
(jGirardi et al.ll2002D in the range 0.7-2 Gyr (Z=0.019). 
We find that the Aa components are consistent with the 
isochrones, but the Ab system components appear over- 
luminous by ~0.5 mag (Figure [5]). However, this must 
be considered a preliminary finding since it is dependent 
on a single measurement of the Aa-Ab magnitude differ- 
ence, and we have not taken into account the (presum- 
ably minor) effects of the slightly different bandpasses 
{K vs. K'). 

3.4. Spectral Classes and v sini 

IStrassmeier fc Fekel ()1990() identified several 
luminosity-sensitive and temperature-sensitive line 
ratios in the 6430-6465 A region. Those critical line 
ratios and the general appearance of the spectrum were 
employed as spectral-type criteria. However, for stars 



that are hotter than about early-G, the line ratios in 
the 6430 A region have little sensitivity to luminosity, 
so only the spectral class of an A or F star can be 
determined. The luminosity class is found by computing 
the absolute visual magnitude with the Hipparcos or 
our orbital parallax and comparing that magnitude to 
evolutionary tracks or a table of canonical values for 
giants and dwarfs. 

The red wavelength spectrum of 88 Tau was com- 
pared with those of slowly rotating Am, F and G dwarfs. 
T hese reference stars, i dentified mostly from the lists 
of lAbt fc MoTTdil (Il995h . iKeenan fc McN^ (fT989h . and 
iFekeil ([l997f ). were observed at KPNO with the same 
telescope, spectrograph, and detector as our spectra 
of 88 Tau. With a computer program developed by 
iHuenemoerder fc BardenI (jl984l ) and. Harden (1985,), var- 
ious combinations of reference-star spectra were rota- 
tionally broadened, shifted in radial velocity, appropri- 
ately weighted, and added together in an attempt to re- 
produce the spectrum of 88 Tau in the 6430 A region. 
lAbt fc Morrelll (|1995f) classified 88 Tau as an Am star 
with spectral classes of A4, A6, and A7 for the calcium, 
hydrogen, and metal lines, respectively. Their classifi- 
cation of HR 3526 was identical to that of 88 Tau, so 
we adopted the spectrum of HR 3526 as the proxy for 
the primary of the 3.57 day binary, which dominates the 
spectrum at blue wavelengths and is still the strongest 
component in our red- wavelength region (Fig. 2). The 
3 subclass difference between the calcium K line type 
and metal line typ e indic ates that this star is a marginal 
or mild Am star (jAbt fc Bidclman 1969). A good fit 
to the fines of the 3.57 day secon dary was produced 
by P rocyon, spectral type F5 IV-V (jJohnson fc Morganl 
Il953| ). Lines of the two components in the 7.89 day bi- 
nary are similar in strength and rotation and were well 
represented by a sp ectrum of 70 Vir, spectral type G4 V 
(jKeenan fc McN eil 198 9i). A fit with HR 483, spectral 
type G1.5 V pCeenan fc McNeill I1989D . used as a proxy 
for the 7.89 day binary pair, was nearly as good. Thus, 
the spectral classes of the four stars are A6m, F5, G2-3:, 
and G2-3:, where the colon indicates that the spectral 
class is more uncertain than usual because of the weak- 
ness of the lines. As shown by their positions in Figure 8, 
all four stars are dwarfs. The abundances of Procyon, 70 
Vir, and HR 483 are essentially solar, indicating that the 
abundances of the components of 88 Tau, except for the 
Am star, are also close to solar. 

The continuum intensity ratio of our best reference 
star combination, fitted to the spectrum of 88 Tau, is 
0.79:0.11:0.05:0.05. If we adopt the continuum intensity 
ratios as the luminosity ratios at 6430 A, we obtain a 
magnitude difference of 2.1 ± 0.3 between the 3.57 day 
pair, and 2.4 ± 0.3 between the astrometric components, 
where the uncertainties are estimated. The 6430 A wave- 
length is about 0.6 of the way between the center of the 
Johnson V and R bandpa sses. 

With the procedure of Fckel' ()1997t ). we determined 
projected rotational velocities for the 4 components of 
88 Tau from 10 KPNO red-wavelength spectra. For each 
spectrum the full-widths at half-maximum of 2 or 3 un- 
blended lines in the 6430 A region were measured and the 
results averaged for each component. The instrumental 
br oaden i ng wa s removed, and the calibration polynomial 
of lFekel (fl99l was used to convert the resulting broad- 
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Fig. 7. — The measured and modeled radial velocities of the 
88 Tau Aa-Ab system, with the motion due to the Aal-Aa2 and 
Abl-Ab2 systems subtracted. 

ening in angstroms into a total line broadening in km s~^ . 
From Fckel (1997, 2003) we assumed a macroturbulence 
of 0.0 for the Am star, 4 km s~^ for the mid-F star, and 
3 km s~^ for the G stars. The resulting v sin i values are 
37 ± 2 and 17 ± 2 km s^^ for the primary and secondary 
of the 3.57 day binary, respectively, and 5 ± 3 km s"-'^ for 
both components of the 7.89 day binary. Our value for 
the Aal component is con sistent with t he de termination 
of vsini= 36 km s'^ bv lRover et"all (|2002[ ). 

To determine whether the rotational velocities of the 
binary components are synchronized, we assumed that 
the rotational axis of each component is parallel to its 
orbital axis. We then computed the equatorial rotational 
velocity for each component and compared it with its ex- 
pected synchronous velocity, computed with its canonical 



radius (jGravl[l993 ). The resulting equatorial rotational 
velocities are 40, 18, 11, and 11 km s~^, for components 
Aal, Aa2, Abl, and Ab2, respectively. The computed 
synchronous rotational velocities are 24, 19, 6.5, and 6.5 
km s~^, respectively. 

Component Aal, the Am star and primary of the 3.57 
day binary, is the only one of the four components that 
does not have a convective atmosphere. It is also the only 
component that is clearly not rotating synchronously: its 
equatorial rotational velocity is 1.67 times faster than 
synchronous. On the other hand its F5 companion, com- 
ponent Aa2 is rotating synchronously. Because of the 
relatively large uncertainties of our v sin i measurements 
for Abl and Ab2, it is possible that the two components 
of the 7.89 day are also synchronously rotating. 

4. CONCLUSION 
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PHASES interferometric astrometry has been used to- 
gether with radial velocity data to measure the orbital 
parameters of the quadruple star system 88 Tau A, and 
in particular to resolve the apparent orbital motion of 
the close Aal-Aa2 and Abl-Ab2 pairs. We have made 
the first determination of the period of the Ab binary 
system and found it to consist of a pair of nearly equal- 
mass G stars. The amplitude of the Abl-Ab2 Center- 
of-Light motion is only ^ 65//as, indicating the level of 
astrometric precision attainable with interferometric as- 
trometry. We are able to resolve the orbital motion of 
all of the components, and hence determine the orbital 
inclinations and component masses with a precision of a 
few percent. Finally, we are able to determine the mutual 
inclinations of the various orbits. 
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